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We show the existence of an exact ground state in certain parameter regimes of one-dimensional 
half-filled extended Hubbard model with site-off-diagonal interactions. In this ground state, the 
bond-located spin correlation exhibits a long-range order. In the case when the spin space is SU(2) 
symmetric, this ground state degenerates with higher spin states including a fully ferromagnetic 
state. We also discuss the relation between the exact bond-ordered ground state and the critical 
bond-spin-density-wave phase. 



The Hubbard model is one of the generic models to de- 
scribe interacting electrons in narrow-band systems |Q. 
The on-site repulsion of this model is due to the matrix el- 
ements of the Coulomb interaction corresponding to the 
on-site Wannier states, and the other matrix elements 
are neglected. The importance of the neglected nearest- 
neighbor exchange interactions was, however, stressed 
and discussed for stabilization of ferromagnetism or a 
dimerized state [|||]. 

Japaridze first discussed using the weak-coupling 
theory the possibility of the bond-spin-density-wave 
(BSDW) ground state connected with the site-off- 
diagonal nature of the pair- hopping term [Q. He intro- 
duced the BSDW order parameters d = {0^,0^, O-): 

<j<j' 

with the Pauli matrices r — (t x , t v , t z ), where the oper- 
ator c\ a (c irT ) creates (annihilates) an electron with spin 
cr (=|, i) at the i-th site. The z-component of these order 
parameters (BSDW- 2: ) Of describes a staggered magne- 
tization with spins located on bonds between adjacent 
sites (see Fig. 0). 

t M M 

3OOOO0 E 

FIG. 1. The BNeel state given by Eq. (|). The enclosed 
pairs of sites (gray ovals) indicate electron-hole dimers. 

First in this Letter, we exactly demonstrate that the 
half-filled one-dimensional (ID) Hubbard model with 
site-off-diagonal interactions and a spin anisotropy has a 
doubly degenerate bond-ordered ground state in an inter- 
mediate coupling regime, using the optimal ground state 
approach |^-^. This bond-ordered state is interpreted 
as a Ncel state of bond-located spins, and the state real- 
izes the fully BSDW long-range order with respect to the 



z-direction. For this reason we call the state the Bond- 
Neel (BNeel) state in the following. In this state there 
are no correlations with respect to both of site-located 
charge and spin density operators except for the nearest 
neighbor. The system possesses, therefore, charge and 
spin gaps between the ground state and excited states. 

Next we consider a spin SU(2) symmetric case, and we 
show that the existence of the phase in which the BSDW 
correlation is dominant by using numerical methods as 
the level-crossing approach |10| , |ll]| . The BSDW correla- 
tion shows a power-law decay at large distances and the 
system has a gapless spin excitation. This BSDW phase 
borders on the BNeel phase on the multicritical point. 
We discuss the relation between the BNeel state and the 
critical BSDW state. 

We consider the ID extended Hubbard model with on- 
site and nearest-neighbor interaction given by 

lEE^I"'"''! C I<7 c ]<y> c i<y' c ka < 

ijkl ijg 1 

where the matrix elements with respect to on-site or 
nearest-neighbor only are nonzero. The Hamiltonian is 
H = hij, where the local Hamiltonian hij is given 

by 

h ij = ^( c L c 3 a + h - c -) H + X ( n i-<r + n ]-<?)] 

+ h u [i n K - h)( n ii ~ 2) + K't - h)( n n - I)] 

cr cr' 

+ iiyy^( c t c + c t \( t + t \ ( 2 ) 

aa' 

Here the summation is taken over pairs of neighbor- 
ing sites and the number operators is defined as 
n ia = c\ a c i(y . The first term in the Hamiltonian (|^) is 
the single-particle hopping (t-term) including the bond- 
charge interaction X = (ii\v aa i \ij). The on-site coupling 
is U = (ii\v a> - a \ii) . Nearest- neighbor density interac- 
tion with opposite and parallel spins are parametrized by 
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Va.-a = V± = {ij\v a - a \ij) and V aa = V\\ = {ij\v aa \ij), 
respectively. Here we introduce spin anisotropy due to 
some factors, such as the crystal structure and spin-orbit 
interaction. The bond-bond interaction is controlled by 
W = (ij\v aa *\ji) = (ii\v aa /\jj). 

Now let us show that the above Hamiltonian can be 
factorized in some parameter space. For this purpose, 
we introduce a local bond operator with spin a defined 
by 

Rij (A, m) = A a\ jrT a ijc + fi s lj<7 s\ JtT , (3) 

where the operators s\j a and aJ J(T denote the symmetric 
and antisymmetric bond-creation operators, respectively: 
4j„ = (4 + a% a = (4 - 4)/V2. The bond 

operator (A, /i) can be rewritten by using the opera- 
tors Bfj = J ia c ]a + h.c. and AT? = n la + n ja : 

= ±±±(1 - + ^(^5 - !)• ( 4 ) 
The local Hamiltonian (^) can be written as follows: 

/iy = 4 ( a , (3) 4 («', /?') + 4 (A, /i) 4 (V, //) 

+ ^K^(n lCT -i)(n^ -i) (5) 

+ (W + X - t) Bfj + X(m + n ) - (t + X), 

a 

where V aa = V\\-W and V aa , = Vj_-U/2 (a == a'). Here 
the parameters of operators Rfj in (|5|) are constrained to 
satisfy the relations: 

2(W + AA') = U + 6W -At, 

2(pfi + lili') = U - 2W + At, (6) 

2(a/3' + A//) = 2(a'/3 + A'/i) = -U + 2W. 

When A, fi > 0, we see from the Eq. (|) that the bond 
operator is positive semidefinite: (i?^ (A,/i)) > 0. Then 

for U = 2V±, W — V\\ = t — X the local operator hij = 
hij — X(m + rij) + 1 + X is also positive semidefinite as 
long as a, a', (3, (3 1 , A, A', fx, li! > 0. 

From these facts, we find that the energy of the global 
Hamiltonian TL is bounded from below in the following 
regime: 

U = 2V±, W = V\i=t-X, 

max{2W - At, -6W + At} < U < 2W. (7) 

We denote the global Hamiltonian H in this regime as 
TL\,. We impose the periodic boundary condition on the 
system and choose the number of lattice sites L to be even 
in the following. The lower bound is given as Ei owcr = 
— (t + X)L + 2XN, where N denotes the total number of 
electrons. To obtain an upper bound of the ground-state 
energy, we consider the following trial wave functions: 



i*5)=n4v n 4-+i,-j°>. ( g ) 

iE{odd} j'G{cvcn} 

where |0) denotes the vacuum. A schematic illustra- 
tion of this state is given in Fig. [j] The state corre- 
sponds to a half-filled state with a density N/L = 1. 
Since Rf 5 (A, fj,) s] |0) = 0, the state sjj-jO) is an eigen- 
state of the local Hamiltonian with eigenvalue zero: 
hij s lja\fy = 0. Rewriting the global Hamiltonian as 
Hb = [hij + X(rii + rij) - (t + X)], we can see 

that the state (||) is also an eigenstate: 

H h \^ ) = -(t-X)L\%). (9) 

This means that the ground-state energy is bounded from 
above at E uppcl = — (t — X)L. At half filling, therefore, 
the upper and lower bound coincide, and the eigenstate 
I^/q) is the ground state in the regime (0). The ground- 
state energy is given as E = — (t — X)L. 
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FIG. 2. Energies in the L = N — 8 system at V = W = t. 
In the region where the singlet excitation (dotted line) is 
lower than the triplet one, the spin excitation has a gap. 
The singlet-triplet level crossing point corresponds to the 
CDW-BSDW phase boundary. 

When U — 2W > t, the system is in the the regime (0), 
and the spin exchange interaction obtains the SU(2) sym- 
metry, Vj| = Vx = V. Then the operators S" ot = J2i i 
a = +, — , z commute with the Hamiltonian. In this case, 
we can deduce that the ground state |^o) an d a fully po- 
larized ferromagnetic (FM) state, |FM) = JT^ cL |0), are 
degenerate. Since \^o) does not vanish, the states 

h h ' ")° = Xl(sL) lk (s^r \n) (io) 

' k 

are also eigenstates with the lowest energy eigenvalue 
and with S* ot = J2k( l k ~ rn k ), (|S? ot | < L/2) as long 

as I 12 " ) CT ^ 0. These states, however, may not 

I mi ni2 ... I 1 1 i j 

be linearly independent each other. The eigenstates of 
S\ ot with Sl ot may be constructed by linear combina- 
tions of the states ( |To| ) with 5f ot . Then these S^ot states 
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(|Stot| — ^/2) may degenerate. In Fig. || we show en- 
ergies in the L = 8 system at W — V — t in Eq. (||) 
using the exact diagonahzation. At the exactly solv- 
able point (U = 2t), all the lowest states with total spin 
I Stot | < L/2 and the density N/ L — 1 degenerate. 

For V± > Vy, i.e., t < 2W < U, the Hamiltonian can 
be expressed as 

« = n h \ u=2W + 7 E K- - !) (4 - !) ( n ) 

with 7 = U/2-W > 0. The ground state is the FM state, 
since |FM) gives the lowest eigenvalue both of 'Hb\ u _ 2W 

and J2(ij)( N lj - 1 )(4 _ !)• Tne system undergoes, 
therefore, a first-order phase transition at this level- 
crossing point 7 = 0. At U = -2W = -2t, X = 0, the 77- 
paring operator with momentum 7T, 77^ = l) J ' c Jl c it> 
commutes with the Hamiltonian. Therefore the 77-paring 
state: \rj) = (f?^)^ 2 ^) is also an exact ground state, 
which agrees with the result by de Boer and Schadschnei- 
der ](|. From the above discussion, we obtain the phase 
diagram of the model (||) as shown in Fig. |^. 
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FIG. 3. The phase diagram of the present model (0) in 
the U-W plain. The parameters of the model are fixed as 
U = 2V± and W = Vj| = t - X . The X-term vanishes on the 
vertical dotted line. The point denotes at V± = Vu,X = 0. 

In order to analyze the properties of two ground states 
|^o) an d |^o) m the regime (|7]), we introduce the matrix 
product representation |7],|| . This is useful for calculating 
the ground state expectation values and correlation func- 
tions. First, we describe these states in terms of products 
of local 2x2 matrices as, 

= 2- L ' 2 Tr gl ® g 2 a ® • • • ® g" L _ x ® g~ L ° , (12) 

where two matrices gf , a =|, |, are given by 

, - ( l 2 >* ^ (13) 

5l { \o)i \-*)< J ■ { 3) 



Here \a)i = cTJO);, |2) 4 = c^cjjO);, and ® denotes 
the usual matrix multiplication of 2 x 2 matrices with 
a tensor product of the matrix elements. Note that 
hij(gf ® 97) = 0- Next, we introduce 4x4 transfer 
matrices: Gf£ = flgt„ 2 fl&' ima , where Z fc = {n k ,m k ) and 
the indices correspond as = 1,2, 3, 4 «-> (nfe,rrife) 

1 1 . (12), (21), (22), respectively. Since the overlap be- 
tween \^q) and |\&o) f° r si ze is evaluated by (^Jl^o) = 
2~ L Tr(G n ) L = 4-2~ L , these two states are orthogonal 
in the limit L — > 00. The ground state is, therefore, dou- 
bly degenerate, and the translational symmetry of the 
system is spontaneously broken except for U = 2W. 

The two-site correlation functions of operators A\ at 
the 1st site and B r at the r-th site can be calculated 
by {A\B r ) = Tr Z(A) G r ~ 2 Z(B) G L ~ r with Z{A) = 
g a \Ag a . Thus, the nearest neighbor correlations are 
obtained as follows, (n.jn i+1 ) = 2, {S*Sf +1 ) = —\, 
and (S^S^ +1 ) = 0. On the other hand, the two-point 
charge and spin correlation functions for \i — j\ > 2 are 
(niUj) — (rii}(nj) = and (S?Sj) = 0, respectively. 
These results indicate that there is a finite energy gap 
between the ground state and the excited states with re- 
spect to site-located charges and spins. 

We can also obtain the bond-bond correlation func- 
tions as (Bf t+1 Bfj +1 ) = (1 ± (-l)l i -J'l)/2, and thus we 
obtain the BSDW-z correlation as (OfO?) = 1. 

This result shows the existence of the BSDW-z long- 
range order, and the expectation values of the BSDW- 
z order parameter is given by (Of) = ±1. Since 

{ s i i+i o-' a ] 3+1 tr} = 0> the expectation value of Of can 
be expressed as (Of) = ((-l) l (M, T - M u )) with M ia = 
s\ i+1 a s i i+1 a . Now we introduce a bond-spin vector op- 
erator given as T l = f E^i+i/^Su+i,,- Then 
we have 

(Ol)=2({-l) i T?), (14) 

The BSDW-z order just corresponds to the Neel order of 
the bond-located spins. Since the nearest-neighbor com- 
mutation relations of the bond-spin operators are differ- 
ent from those of the ordinary spin operators, the bond- 
spin and the spin operators are not exactly equivalent. 

Next we discuss the relation between the BNeel and 
the BSDW states. As was discussed above, the BNeel 
state has both charge and spin gaps, whereas the BSDW 
state has gapped charge and gapless spin excitations [Q. 
Let us consider the Hamiltonian with an SU(2) symmetry 
in spin space given at X = 0, V± = V\\ = V in Eq. (^). In 
this Hamiltonian, j^o) gives an exact ground state when 
U = 2V — 2W — 2t. There is large degeneracy in the 
ground state and the correlation function of the BNeel 
state are not reliable. 

We show in Fig. |J the phase diagram at W/t = 1 in 
the U-V plain, obtained by the exact diagonahzation 
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of the L = 12 system. There appear BSDW, charge- 
density- wave (CDW), and FM phases. The CDW state 
has gaps both in charge and spin excitations, and shows a 
site-long-range order. The spin-gap transition occurs be- 
tween the CDW-BSDW boundary. This transition point 
is obtained by the singlet-triplet level crossing in the 
excited states, which is justified by the conformal field 
theory @. On the other hand, the CDW-FM and the 
BSDW-FM transitions are of the first order, and their 
boundaries are determined by level crossing in the ground 
state O]. Finite size effects in this analysis are small 
enough to get reliable results. The exact ground state 
appears just on the BSDW-FM boundary within the pre- 
cision of the exact diagonalization. 



In the case when SU(2) symmetry exists in the spin 
space, the ground state degenerates with higher spin 
states including the fully FM state. This BNeel phase 
adjoins the critical BSDW phase on the multicritical 
point. We expect that the relation between the BNeel 
and the BSDW states is analogous to the one between 
the Ising limit and the Heisenberg point in the spin-1/2 
antiferromagnetic XXZ chain with respect to the low- 
energy spin excitations. Note, in the model (^), that the 
BNeel state appears when the spin anisotropy is XY-likc, 
V±/W = V±/V\\ < 1, contrary to the Neel state of the 
XXZ model. 
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FIG. 4. The phase diagram of the model in the SU(2) sym- 
metric case with X = 0, W = t, obtained by the numerical 
data of the L = 12 system. The wave function of Eq. (^) gives 
the exact ground state at U = 2V = 2t, which corresponds to 
the point indicated in Fig. |j| 

As was shown in Fig. ||, at the exactly solvable point 
(U = 2t), all the lowest states with the total spin Stot 
(Stat < L/2) degenerate. From the above numerical 
and exact results, we conjecture that this total spin-S , to t 
states can be constructed by means of linear combina- 
tions of the states ( fio| ) and the singlet state in the con- 
structed states corresponds to the BSDW state. 

Finally, we comment on the relation between the 
present result and other works. Exact dimerized ground 
states in electron systems is discussed recently by 
Dmitriev et al. |jr2f . The dimer consists of up- and down- 
spin electrons, whereas our "dimer" given in Eq. (|^) of 
an electron and a hole for each spin. Their results, there- 
fore, do not contain a state with the staggered magneti- 
zation on the bonds. The ground state we have discussed 
is rather similar to the one which appears in a spin-1/2 
two-leg ladder model H or in an orbitally degenerate spin 
chain ]j| with the Jordan- Wigner transformation. 

In summary, we have shown that a BNeel state with 
the BSDW order is the exact ground state in wide pa- 
rameter space of the generalized Hubbard chain. In 
this state, both charge and spin excitations have gaps. 
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